Computational  Modeling  of  Sense 
and  Destroy  Armor  (SAD ARM) 
Submunition  Separation/Collision 

by  Jubaraj  Sahu,  Karen  R.  Heavey, 
and  Charles  J.  Nietubicz 


ARL-TR-1378 


June  1997 


19970630  125 


Approved  for  public  rdease;  distribution  is  unlimited. 


D2TC 


The  findings  in  this  report  are  not  to  be  construed  as  an  ofi&dal 
Departmrait  of  the  Army  position  unless  so  deagnated  by  other 
auliiorized  documents. 

Citation  of  manu&cturiar’s  a:  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  thereof. 

Destroy  ftis  rqxjrt  when  it  is  no  longa:  need.  Do  not  return  it  to 
the  originator. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5066 


ARL-TR-1378  June  1997 


Computational  Modeling  of  Sense  and 
Destroy  Armor  (SAD ARM)  Submunition 
Separation/  Collision 


Jubaraj  Sahu,  Karen  R.  Heavey,  Charles  J.  Nietubicz 

Weapons  and  Materials  Research  Directorate,  ARL 


Approved  for  public  release;  distribution  is  unlimited. 


Abstract 


Sense  and  Destroy  Armor  (SADARM)  is  a  high-priority  Army  program  that  experienced  an 
initial  design  flaw,  leading  to  submunition  collisions.  Computational  fluid  dynamics  (CFD), 
including  Chimera  overset  technology,  has  recently  been  applied  to  solve  this  problem.  The 
CFD  modeling  has  provided  significant  physical  insight  into  the  complex  unsteady  flow  field 
associated  with  the  SADARM  submunitions  during  the  separation/ejection  process.  CFD 
simulations  have  predicted  that,  with  the  addition  of  the  fins  on  the  trailing  submxmition, 
separation  will  occur,  and  collisions  will  be  avoided. 
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1.  INTRODUCTION 


An  important  parameter  in  the  design  of  shell  and  bodies  flying  in  relative  motion  to  each 
other  is  the  total  aerodynamic  drag.  The  base  drag  constitutes  a  large  part  of  the  total 
aerodynamic  drag,  and  accurate  prediction  of  the  base  region  flow  field  is  necessary.  The  ability 
to  compute  the  base-region  flow  field  for  projectile  configurations  using  Navier-Stokes 
computational  techniques  has  been  developed  over  the  past  several  years  [1,2].  Recently, 
improved  numerical  predictions  have  been  obtained  using  a  more  advanced  zonal  upwind  flux- 
split  algorithm  [3].  This  zonal  scheme  preserves  the  base  comer,  which  allows  an  accurate 
modeling  of  the  base-region  flow.  Previous  computational  studies  have  been  completed  showing 
the  aerodynamic  effect  for  a  variety  of  base  geometries.  These  calculations,  however,  were 
performed  on  stand-alone  projectile  configurations  and  represent  a  single-body  problem. 
Recently,  a  multibody  problem  that  involves  other  bodies  flying  in  the  wake  of  a  parent  projectile 
has  required  computational  analyses  [4].  This  is  due  in  part  to  the  difficulty  in  finding  good 
experimental  and/or  analytical  data  for  such  problems.  The  particular  problem  here  was  to 
determine  the  aerodynamic  effect  of  small  cylindrical  segments  being  ejected  into  the  wake  of  a 
parent  projectile.  The  complexity  and  uniqueness  of  this  problem  resulted  from  the  trailing 
segments  being  in  relative  motion  to  each  other,  embedded  in  a  nonuniform  wake  flow,  and 
requiring  a  time-dependent  solution.  These  calculations  were  performed  at  a  supersonic  speed 
for  the  dynamic  case,  which  involved  time-accurate  numerical  computations,  and  the  computed 
results  were  foimd  to  be  in  good  agreement  with  the  test  data. 

The  present  problem  is  another  multibody  problem  that  involves  separation  of  two 
submunitions  at  a  low  transonic  speed.  The  Sense  and  Destroy  Armor  (SAD ARM)  program  is  a 
high-priority  Army  program  that  had  experienced  an  initial  design  flaw  leading  to  submunition 
collisions.  Figure  1  shows  the  components  for  the  initial  system.  S^aration  of  the  submunitions 
depends  on  the  initial  ejection  velocity.  For  small  initial  ejection  or  separation  velocity,  the 
submunitions  can  in  fact  s^arate  initially,  but  later  in  the  flight  are  pulled  back  together  to 
collide.  Figure  2  shows  the  collision  of  two  submunitions  in  a  flight  test  A  time-accurate 
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Figure 


Figure  2. 


numerical  solution  of  this  problem  is  the  subject  matter  of  this  research  report,  intended  to 
provide  more  basic  flow-field  information  and  understanding  to  avoid  such  collisions. 

The  time  accurate  numerical  simulation  of  the  multiple  aerodynamic  bodies  in  relative 
motion  has  been  obtained  using  the  CSiimera  [5]  approach.  This  technique  has  been  used  to 
compute  inviscid  and  viscous  flows  about  complex  configurations  [6, 7],  and  has  been 
demonstrated  for  unsteady  viscous-flow  problems  with  bodies  in  relative  motion  [8].  The 
Chimera  approach  is  a  domain  decomposition  method  that  uses  overset,  body-conforming  grids 
and  grew  out  of  the  necessity  to  computationally  model  geometrically  complex  configurations. 
The  originally  developed  Navier-Stokes  numerical  capability,  zonal  F3D  [3],  was  extended  by 
the  authors  to  include  the  details  of  the  Chimera  procedure.  This  work  couples  the  solution  of 
the  Navier-Stokes  equations,  which  govern  fluid  motion,  with  the  solution  to  the  6  degree-of- 
freedom  (DOF)  equations  of  motion.  The  coupling  of  the  fluid-dynamic  solution  and  rigid-body 
motion  eliminates  the  need  for  simplifying  assumptions  and  allows  more  accurate  physically 
based  simulations. 

2.  GOVERNING  EQUATIONS  AND  SOLUTION  TECHNIQUE 

The  complete  set  of  time-dependent,  Reynolds-averaged,  thin-layer  Navier-Stokes  equations 
is  solved  numerically  to  obtain  a  solution  to  this  problem.  The  numerical  technique  used  is  an 
implicit,  finite-difference  scheme.  Time-accurate  calculations  are  made  to  numerically  simulate 
the  separation  and  collision  of  two  submunitions. 

2.1  Governing  Equations.  The  complete  set  of  three-dimensional  (3-D),  time-dependent, 
generalized  geometry,  Reynolds-averaged,  thin-layer  Navier-Stokes  equations  for  general  spatial 
coordinates  Ti,  and  C  can  be  written  as  follows  [9]: 

a^q  +  +  3^6  +  djH  =  Re-^a^S,  (1) 

where 

$  =  5(x,  y,z,  t)  -  longitudinal  coordinate; 

q  -  T|(x,  y,  z,  t)  -  circumferential  coordinate; 

C  =  C(x,y,  z,  t)  -  nearly  normal  coordinate; 

T  -  t  -  time 
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In  Equation  1,  the  thin-layer  approximation  is  used,  and  the  viscous  terms  involving  velocity 
gradients  in  both  the  longitudinal  and  circumferential  directions  are  neglected.  The  viscous 
terms  are  retained  in  the  normal  direction,  C,  for  the  projectile  and  segments,  and  are  collected 
into  the  vector  S .  hi  the  wake  or  the  base  region,  similar  viscous  terms  are  also  added  in  the 
streamwise  direction,  For  this  computation,  the  diffusion  coefficients  p  and  k  contain 
molecular  and  turbulent  parts.  The  turbulent  contributions  are  supplied  through  an  algebraic 
eddy  viscosity  turbulence  model  developed  by  Baldwin  and  Lomax  [10]. 

The  velocities  in  the  t),  and  C  coordinate  directions  can  be  written  as 

U  =  ^t  +  u5x  +  v5y  +  w^2 
V  =  llt  +  UTl^  +  VTly  +  Wll^ 

W  =  Ct  +  uCx  +  vCy  +  wC^, 

which  represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  retained  as  the  dependent  variables  and  are 
nondimensionalized  with  respect  to  a„  (the  ffee-stream  speed  of  sound).  The  local  pressure  is 
determined  using  the  relation 

P  =  (Y  “  1)  [e  -  0.5p  (u^  +  v^  +  w^)] ,  (4) 

where  y  is  the  ratio  of  specific  heats.  Density,  p,  is  referenced  to  p„  and  the  total  energy,  e,  to  p„ 
aj.  The  transport  coefficients  are  also  nondimensionalized  with  respect  to  the  corresponding 

A 

free-stream  variables.  Thus  the  Prandtl  number  that  appears  in  S  is  defined  as  Pr  = 

2.2  Numerical  Technique.  The  implicit,  approximately  factored  scheme  for  the  thin-layer 
Navier-Stokes  equations  using  central  differencing  in  the  tj  and  C  directions  and  upwinding  in  ^ 
is  written  in  the  following  form  [11], 
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I  +  i^h6^(A+)“  +  -  i^Dj  1^ 

X  [l  +  i^h5|(A-)"  +  i,h6,6"  - 

=  -i,At{6;‘[(Fi"  -  ^:]  +  5|[(^-)"  -  ^.]  *  6^6”  -  &.) 

+  6j(h«  -  H.)  -  Re->8((§”  -§.)}-  ijD^Q"  -  Q.). 


(5) 


where  h  =  At  and  the  free-stream  base  solution  is  used.  Here,  6  is  typically  a  three-point  second- 
order  accurate  central  difference  operator,  6  is  a  midpoint  operator  used  with  the  viscous  terms, 
and  the  operators  6^*’  and  6^^  are  backward  and  forward  three-point  difference  operators.  The 
flux  F  has  been  eigensplit,  and  the  matrices  A,  B,  C,  and  M  result  from  local  linearization  of 
the  fluxes  about  the  previous  time  level.  Here,  J  denotes  the  Jacobian  of  the  coordinate 
transformation.  Dissipation  operators  and  Dj  are  used  in  the  central  space  differencing 
directions.  The  smoothing  terms  used  in  the  present  study  are  of  the  form: 


DJ,  =(At)J-‘ 


ej8p(B)p8 


P(B)  ^3 
1+p 


Djl,  =  (At)I-'[e25p(B)p8  ♦  2,5e,8p(B)S]|,J. 


where 


P  = 


_ yp| 

1  (1  +  6^)P 


and  where  p(B)  is  the  true  spectral  radius  of  B.  The  idea  here  is  that  the  fourth  difference  will  be 
tuned  down  near  shocks  (e.g.,  as  P  gets  large,  the  weight  on  the  fourth  difference  drops  down 
while  the  second  difference  tunes  up). 
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2.3  rhimera  Composite  Grid  Scheme.  The  Chimera  overset  grid  scheme  is  a  domain 
decomposition  approach  where  a  configuration  is  meshed  using  a  collection  of  overset  grids.  It 
allows  each  component  of  the  configuration  to  be  gridded  separately  and  overset  into  a  main  grid. 
Overset  grids  are  not  required  to  join  in  any  special  way.  Usually  there  is  a  major  grid  that 
covers  the  entire  domain  or  a  grid  generated  about  a  dominant  body.  Minor  grids  are  generated 
about  the  rest  of  the  other  bodies.  Because  each  component  grid  is  generated  independently, 
portions  of  one  grid  may  be  found  to  lie  within  the  solid  boundary  contained  within  another  grid. 
Such  points  lie  outside  the  computational  domain  and  are  excluded  from  the  solution  process. 

In  the  submunition  separation  problem,  the  first  (leading)  submunition  is  a  major  grid  and  the 
second  (trailing)  submunition  grid  is  a  minor  grid.  The  minor  grid  is  completely  overlapped  by 
the  major  grid;  thus,  its  outer  boundary  can  obtain  information  by  interpolation  from  the  major 
grid.  Similar  data  transfer  or  communication  is  needed  from  the  minor  grid  to  the  major  grid. 
However,  a  natural  outer  boundary  that  overlaps  the  two  submunition  grids  does  not  exist.  The 
Chimera  technique  creates  an  artificial  boundary  (also  known  as  a  hole  boundary)  between  grids 
that  provides  the  required  path  for  information  transfer  from  the  first  submunition  to  the  second. 
The  resulting  hole  region  is  excluded  from  the  flow-field  solution  in  the  projectile  grid.  Equation 
5  has  been  modified  for  Chimera  overset  grids  by  the  introduction  of  the  flag  ibto  achieve  just 
that.  This  ij,  array  accommodates  the  possibility  of  having  arbitrary  holes  in  the  grid.  The  i^  array 
is  defined  such  that  ib  =  1  at  normal  grid  points  and  ib  =  0  at  hole  points.  Thus,  when  ib  =  1, 
equation  5  becomes  the  standard  scheme.  But,  when  ib  =  0,  the  algorithm  reduces  to 
AQ”  =  0  or  ^  ^  =  Q“ ,  leaving  Qunchanged  at  hole  points.  The  set  of  grid  points  that  form 
the  border  between  the  hole  points  and  the  normal  field  points  are  called  intergrid  boundary 
points.  These  points  are  updated  by  interpolating  the  solution  from  the  overset  grid  fliat  created 
the  hole.  Values  of  the  ib  array  and  the  interpolation  coefficients  needed  for  this  update  are 
provided  by  a  separate  algorithm  [6]. 

In  die  present  study,  which  involves  multiple  bodies  in  relative  motion,  the  location  of  the 
holes  and  the  intergrid  boimdary  points  are  time  dependent.  Accordingly,  the  ib  array  and  the 
interpolation  coefficients  are  functions  of  time.  This  procedure  of  unsteady  Chimera 
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decomposition  has  been  successfully  demonstrated  by  Meakin  [8].  The  method  depends  on  three 
functions:  domain  connectivity,  aerodynamics,  and  body  dynamics.  The  aerodynamics  code 
depends  on  the  domain  connectivity  code  to  supply  hole  and  interpolation  information.  The 
domain  connectivity  code,  in  turn,  depends  on  the  body  dynamics  code  to  supply  the  location  and 
orientation  of  the  moving  bodies  relative  to  the  primary  body.  Finally,  the  body  dynamics  code 
depends  on  the  aerodynamics  code  to  provide  the  aerodynamic  forces  and  moments  acting  on  the 
moving  bodies. 

The  Chimera  procedure  reduces  a  complex  multibody  problem  into  a  number  of  simpler 
subproblems.  For  moving-body  problems,  all  grids  are  allowed  to  move  with  6  DOF  relative  to 
an  inertial  reference  frame.  Accordingly,  bodies  can  move  with  respect  to  others  without  the 
necessity  of  generating  new  grids.  With  this  composite  overset  grid  approach,  it  is  thus  possible 
to  determine  the  unsteady  relative  motion  of  the  segments  and  associated  aerodynamic  forces 
without  the  need  for  costly  regridding.  This  also  eliminates  potential  accuracy  problems  due  to 
severe  grid  stretching  used  by  many  other  techniques.  The  solution  procedure  is  to  compute  the 
flow  field  at  each  time  step;  integrate  the  pressure  and  viscous  forces  for  the  trailing  segments  to 
obtain  the  drag  force;  use  the  predicted  drag  in  a  coupled  6-DOF  program  to  compute  the  new 
relative  position  of  the  trailing  segment.  Because  this  is  an  axisymmetric  calculation  (0°  angle  of 
attack),  only  the  drag  coefficient  is  used  to  determine  the  motion  of  the  segments.  At  the  next 
time  step,  the  solution  procedure  is  repeated  for  the  new  position  with  the  domain  decomposition 
providing  all  the  required  interpolation  information.  Computations  are  performed  on  each  grid 
separately.  These  grids  use  the  available  core  memory  one  grid  at  a  time.  The  remaining  grids 
are  stored  on  an  external  disk  storage  device  such  as  the  solid-state  disk  (SSD)  device  of  the  Cray 
C-90  or  Y-MP  computer. 

2.4  PomaiTi  Connectivity  Function.  A  major  part  of  the  Chimera  overset  grid  approach  is 
the  information  transfer  from  one  grid  into  another  by  means  of  the  intergrid  boundary  points. 
Again,  diese  points  consist  of  a  set  of  points  that  define  the  hole  boundaries  and  outer  boundaries 
of  the  minor  grids.  These  points  depend  on  die  solutions  in  the  overlapping  regions.  In  the 
present  work,  the  Domain  Connectivity  Function  in  3-D  (DCF3D)  Code  [8]  has  been  used  to 
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establish  the  linkages  between  the  various  grids  that  are  required  by  the  flow  solver  or 
aerodynamics  code  described  earlier.  These  include  the  determination  of  the  interpolation 
coefficients  and  the  setting  up  of  Chimera  logic  for  bodies  making  holes  in  overlapping  grids. 

For  unsteady  moving-grid  cases,  this  code  must  be  executed  at  each  time  iteration.  To  minimize 
the  computation  time,  this  code  uses  the  knowledge  of  hole  and  interpolated  boundary  points  at 
time  level  n  to  limit  its  search  regions  for  finding  their  corresponding  locations  at  time  level  n+1. 

hi  general,  each  component  grid  in  an  overset  grid  system  represents  a  curvilinear  system  of 
points.  However,  the  positions  of  all  points  in  all  the  grids  are  defined  relative  to  an  inertial 
system  of  reference.  To  provide  domain  connectivity,  inverse  mappings  are  used  that  allow  easy 
conversion  from  the  x,  y,  z  inertial  system  to  C»  ti  computational  space.  For  moving-body 
problems,  these  maps  for  component  grids  are  created  only  once.  Identification  of  the  intergrid 
boundary  points  that  correspond  to  the  outer  boundaries  of  the  minor  grids  is  done  simply  by 
specifying  appropriate  ranges  of  coordinate  indices.  The  rest  of  the  intergrid  boundary  points 
that  result  from  holes  created  by  a  body  in  overset  grids  is  a  little  more  difficult  to  identify.  A 
collection  of  analytical  shapes  such  as  cones,  cylinders,  and  boxes  are  used  to  cut  holes  in  tiiis 
method. 

2.5  Boundary  Conditions.  For  simplicity,  most  of  the  boundary  conditions  have  been 
imposed  explicitly  [3].  An  adiabatic  wall  boundary  condition  is  used  on  die  body  surface,  and 
the  no-slip  boundary  condition  is  used  at  the  wall.  The  pressure  at  the  wall  is  calculated  by 
solving  a  combined  momentum  equation.  Free-stream  boundary  conditions  are  used  at  the 
inflow  boundary  as  well  as  at  the  outer  boundary.  A  symmetry  boundary  condition  is  imposed  at 
the  circumferential  edges  of  the  grid,  while  a  simple  extrapolation  is  used  at  the  downstream 
boundary.  A  combination  of  symmetry  and  extrapolation  boundary  conditions  is  used  at  the 
center  line  (axis).  Since  the  free-stream  flow  is  supersonic,  a  nonreflection  boundary  condition  is 
used  at  the  outer  boundary.  Similar  boundary  conditions  are  used  for  the  segments. 
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3.  MODEL  GEOMETRY  AND  COMPUTATIONAL  GRID 


Each  submunition  is  a  right  circular  cylinder  of  length-to-diameter  ratio  1.38.  Figure  3  shows 
an  expanded  view  of  a  composite  computational  grid  for  these  submunitions.  The  complete  grid 
consists  of  approximately  102,660  grid  points  and  is  split  into  8  sections.  The  grid  in  the  wake 
region  of  each  submunition  consists  of  a  large  number  of  pomts  in  both  the  streamwise  and 
normal  directions.  The  surface  points  for  each  region  (body  and  wake)  are  selected  using  an 
interactive  design  program.  Each  grid  section  was  obtained  separately  and  then  appended  to 
provide  the  full  grid.  The  grid  for  the  body  region  as  well  as  the  wake  region  was  obtained 
algebraically.  A  body-conforming  grid  was  obtained  for  each  submunition  and  then  overset  to 
form  the  composite  mesh  shown  in  Figure  3.  This  case  corresponds  to  the  matted  position  when 
both  submunitions  are  merged  together  (flying  as  one  body).  The  major  grid  corresponds  to  the 
grid  for  the  first  submunition,  which  is  easily  generated  independently  of  the  minor  grid  (grid  for 
the  second  submunition).  The  Chimera  technique,  as  stated  earlier,  allows  individual  grids  to  be 
generated  with  any  grid  topology,  thus  making  the  grid  generation  process  easier.  For 
moving-  body  problems,  both  major  and  minor  grids  shown  in  Figure  3  in  the  matted  position 
move  with  the  submunitions  as  they  separate  from  each  other.  Figure  4  is  a  grid  at  another 
instant  in  time  when  the  minor  grid  has  moved  downstream  into  the  wake  of  the  major  grid. 
Again,  there  is  no  need  to  generate  new  grids  for  the  submunitions  (minor  and  major)  during  the 
dynamic  process. 

4.  RESULTS 

Time-accurate  calculations  have  been  performed  to  numerically  simulate  the  separation  of 
two  submunitions,  both  being  in  relative  motion.  All  computations  were  run  at  M.  =  0*80  and 
a  =  0°,  and  atmospheric  flight  conditions  were  used.  These  calculations  required  15M  words  of 
memory,  and  each  case  used  about  20  hours  of  computer  time  on  the  Cray  C-90  supercomputer. 
In  order  to  numerically  determine  the  submunition  separation/collision  process,  computations 
were  performed  for  several  different  initial  separation  velocities.  Results  are  shown  first  for  a 
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case  where  the  initial  separation  velocity  is  5  ft/s.  Figure  5  shows  the  Mach  number  contours  for 
the  submunitions  at  three  points  in  time:  9, 19.5,  and  33.8  ms  after  separation  began.  As  seen  in 
this  figure,  the  flow  field  is  unsteady  and  the  second  submunition  is  completely  submerged  in  the 
wake  of  the  first  submunition.  The  second  submunition  initially  separates  with  the  separation 
distance  having  increased  in  the  time  interval  from  9.0  ms  to  19.5  ms.  Soon  afterward,  both 
submunitions  pull  back  toward  each  other  with  the  separation  smaller  at  33.3  ms  than  at  9.0  ms. 
The  time-dependent  numerical  solution  simulates  the  motion  of  both  submunitions,  and  the 
aerodynamic  drag  changes  continually  (see  Figure  6).  Surface  pressures  including  the  base 
pressures  and  the  viscous  stresses  are  known  from  the  computed  flow  field  and  are  integrated  to 
give  the  aerodynamic  drag  for  both  cylindrical  submunitions.  The  drag  coefficient  for  the 
submunitions  is  shown  in  Figure  6  as  a  function  of  the  separation  distance.  This  drag  coefficient 
is  based  on  the  cross-sectional  area  of  the  submunition  and  the  free-stream  dynamic  pressure. 

As  mentioned  earlier,  the  second  submunition  is  submerged  in  the  subsonic  wake  of  the  first 
submunition.  The  pressure  behind  the  second  submunition  is  higher  than  the  pressure  ahead  of  it 
and,  therefore,  results  in  negative  drag.  The  fluctuating  drag  for  both  submunitions  indicates  that 
the  flow  field  is  unsteady.  The  drag  for  the  first  submunition  is  positive,  which  tends  to  slow  its 
motion.  The  drag  for  the  second  submunition  is  negative,  which  leads  to  it  being  pulled  back 
towards  the  first  submunition.  As  seen  in  this  figure,  the  separation  distance  increases  to  about 
0.8  inches  and  then  decreases.  When  both  submunitions  come  back  close  to  each  other,  the  drags 
for  these  submunitions  are  found  to  be  different  from  what  they  were  during  the  initial 
separation.  This  is  due  to  the  fact  that  the  flow  field  experienced  by  the  submunitions  during 
these  processes  are  different.  The  drag  contributions  for  each  of  the  submunition  surfaces  (front, 
back,  and  top)  are  shown  in  the  next  two  figures.  In  Figure  7,  both  the  front  and  the  back 
surfaces  of  the  first  submunition  are  contributors,  although  the  front  surface  is  shown  to  be  the 
significant  contributor  to  the  drag.  Figure  8  shows  the  drag  from  the  second  submunition  is 
mainly  from  the  front  surface.  For  both  submunitions,  the  top  surfaces  have  little  impact  on  the 
total  drag.  Figure  9  shows  the  separation  distance  as  a  function  of  time  for  the  initial  separation 
velocity  of  5  ft/s.  As  seen  here,  the  separation  distance  increases  with  time  until  a  peak  is 
reached  at  about  20  ms.  The  separation  distance  is  clearly  seen  to  decrease  from  that  point  on, 
and  the  two  submunitions  come  back  together  and  collide  at  38  ms. 
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Figure  6.  Drag  coefficient.  M  =  0.80.  a  =  0.  V.=  5  ft/s  . 
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Figure?.  Drag  coefficient.  M  =  .80.  a  =  0.  V,=  5  ft/s  ffirst  submunition). 
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Figure  8.  Drag  coefficient.  M  =  0.80.  oc  =  0.  V.  =  5  ft/s  (second  submunition). 
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Figure  9.  Separation  distance  vs.  time.  M  =  0.80.  a  =  0.  V.=  5  ft/s. 
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TniriflUy,  the  calculations  for  the  above  case  were  done  using  a  time  interval  of  0.01,  using 
about  50  hours  of  central  processing  unit  (CPU)  time.  A  second  computation  was  performed 
using  a  value  of  0.03  for  the  time  interval.  This  calculation  used  approximately  20  CPU  hours. 
Figure  10  shows  that  the  drag  histories  for  both  calculations  are  similar  and  reach  the  same 
separation  distance  before  being  pulled  back  towards  each  other.  Similarly,  Figure  1 1  shows  the 
peak  separation  distance  for  both  submunitions  to  be  identical,  and  the  time  of  collision  is  very 
similar.  Based  on  these  findings,  the  remainder  of  the  calculations  were  done  using  the  larger 
time  interval  of  0.03,  providing  a  savings  of  over  50%  of  CPU  time. 

A  modification  to  the  projectile  design  resulted  in  an  increase  in  the  separation  velocity. 
Consequently,  a  new  set  of  calculations  were  performed,  this  time  with  an  initial  separation 
velocity  of  15  ft/s.  Figure  12  shows  the  Mach  number  contours  for  the  submunitions  at  three 
points  in  time:  9, 41,  and  69  ms  after  the  onset  of  separation.  The  separation/collision  process 
follows  the  same  trend  as  demonstrated  by  the  5-ft/s  case,  with  the  submunitions  initially 
separating.  Again,  as  shown  in  Figure  13,  the  drag  for  the  first  submunition  is  positive,  slowing 
its  motion  while  the  negative  drag  of  the  second  submunition  pulls  it  back  towards  the  first 
submunition.  After  approximately  40  ms,  the  peak  separation  distance  of  4  inches  is  reached  and 
then  the  second  submunition  is  pulled  back  toward  the  first  submunition.  Figure  14  clearly 
shows  that  collision  occurs  at  80  ms. 

To  avoid  the  submunition  collision,  it  was  thus  necessary  to  either  increase  the  initial 
separation  velocity  or  incorporate  devices  such  as  fins  to  the  second  submunition  to  provide 
added  drag.  Because  of  system  constraints,  the  initial  separation  velocity  could  not  exceed 
15  ft/s.  Thus,  the  second  option  was  selected  and  studied  both  numerically  and  experimentally. 
Figure  15  shows  the  improved  projectile  design.  Fins  were  added  to  the  second  submunition  as 
well  as  to  the  base  plug,  thus  providing  extra  drag.  A  computational  model  of  a  finned 
submunition  requires  a  3-D  calculation  and  large  amoxmts  of  computer  time  and  resources. 
Because  of  time  constraints,  the  existing  computer  model  was  modified  by  adding  extra  drag  to 
simulate  fin-induced  drag.  The  amount  of  added  drag,  as  well  as  the  time  in  which  the  drag 
devices  deployed,  was  obtained  from  an  experiment  conducted  on  the  submunition  alone.  The 
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Figure  10.  Total  drag  for  two  .submunitions.  M  =  0.80.  k  =  0.  V-=  5  ft/s.  Dt  =  0.01.  and  0.03. 
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Figure  11.  Separation  distance  vs.  time.  M  =  0.80.  k  =  0.  V.=  5  ft/s.  Dt  =  0.01 .  and  0.03. 
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Figure  13.  Drag  coefficient.  M  =  0.80.  «  =  0.  V.  =  15  ft/s. 
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Figure  14.  Separation  distance  vs.  time.  M  =  0.80.  a  =  0.  V.  =  15  ft/s. 


19 


Figure  15.  Improved  design  for  SAD  ARM  projectile. 


Chimera  composite  overset  grid  approach  was  then  used  to  numerically  model  this  experimental, 
time-dependent  separation  process.  Figure  16  shows  the  Mach  number  contours  for  the 
submunitions  at  three  points  in  time  (9, 78,  and  1 17  ms)  after  separation  began.  As  seen  in  this 
figure,  the  second  submunition  is  completely  submerged  in  the  wake  of  the  first  submunition  at  9 
ms  and  begins  to  separate  out  from  this  wake  at  a  later  time.  The  added  drag  takes  effect  at  about 
35  ms  of  flight.  The  second  submunition  is  initially  seen  to  separate  and  then  continues  to 
separate  (does  not  come  back).  The  drag  coefficient  for  the  submunitions  is  shown  in  Figure  17 
as  a  function  of  the  separation  distance.  The  drag  for  the  first  submunition  is  again  positive, 
whereas  initially  before  the  drag  devices  are  employed,  the  second  submunition  experiences 
negative  drag  up  to  a  separation  distance  of  about  8  inches.  The  drag  for  the  first  submunition 
decreases,  and  the  drag  for  the  second  submunition  increases  from  that  time  on  until  a  separation 
distance  of  20  inches  is  reached.  Although  drag  fluctuates  somewhat,  the  mean  stays  at  the  same 
value.  Figure  18  shows  the  separation  distance  as  a  function  of  time  for  this  case.  As  seen  here, 
the  separation  distance  increases  with  time  up  to  40  ms.  It  is  seen  to  level  out  around  40  ms,  and 
then  begins  to  increase  again.  The  separation  distance  is  clearly  seen  to  increase  from  that  point 
on,  and  the  two  submunitions  never  come  back  together.  Subsequently,  a  flight  test  was 
conducted.  The  computed  results  are  compared  in  this  figure  with  the  available  experimental  data 
[12]  shown  in  shaded  circles  and  are  found  to  be  in  good  agreement. 
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Figure  17.  Drag  coefficient.  M  =  0.80.  o'  =  0.  V,  =  15  ft/s  with  added  drag. 
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Figure  18.  Separation  distance  vs.  time.  M  =  0.80.  a  =  0.  V.  =  15  ft/s  with  added  drag. 
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5.  CONCLUDING  REMARKS 


A  computational  study  has  been  undertaken  to  compute  the  time-dependent  separation  of  two 
cylindrical  submunitions  in  motion.  Row-field  computations  have  been  performed  at  the  low 
transonic  Mach  number  M»  =  0*80  and  angle  of  attack  a  =  0.0°  using  an  unsteady,  zonal  F3D 
Navier-Stokes  code  and  the  Chimera  composite  grid  discretization  technique.  The  computed 
results  show  the  qualitative  features  of  the  wake  flow  field  for  both  submunitions  for  three 
different  cases.  The  predicted  flow  field  over  the  submunitions  was  found  to  undergo  significant 
changes  as  the  submunitions  separated  from  each  other.  For  initial  separation  velocities  of  both  5 
ft/s  and  15  ft/s,  the  second  submunition  experiences  negative  drag  for  the  dynamic  case  at  all 
times.  In  both  cases,  the  separation  distance  increases  with  time  initially  until  a  peak  is  reached 
and  then  decreases  until  the  two  submunitions  come  back  and  collide.  Time-dependent 
computations  were  also  performed  to  simulate  another  case  where  drag  devices  were  employed 
after  35  ms  of  flight,  with  an  initial  velocity  of  15  ft/s.  Comparison  of  the  computed  results  for 
this  case  have  been  made  with  available  experimental  results,  and  computed  separation  distances 
were  found  to  be  in  good  agreement  with  the  experimental  results.  Both  numerical  and 
experimental  results  show  continued  separation  of  the  submunitions  and  thus  no  collisions  for 
this  case. 

This  work  represents  a  major  increase  in  capability  for  determining  the  aerodynamics  of 
multibody  configurations  and  indeed  has  been  successfully  applied  to  the  SADARM  program.  It 
has  resulted  in  savings  of  time  and  dollars  and  led  to  the  successful  March  1995  Defense 
Acquisition  Board  action  that  resulted  in  low-rate  production  approval  for  SADARM.  The 
complex  physics  and  fluid  dynamics  structure  of  the  aerodynamic  interference  for  the  SADARM 
multibody  problem  has,  for  the  first  time,  been  uniquely  identified  and  could  not  have  been 
accomplished  without  the  computational  fluid  dynamics  (CFD)  techniques  developed  in  ftiis 
research  effort.  This  developed  technology  is  additionally  applicable  to  other  Army  problems 
such  as  time-dependent  sabot  separation  for  kinetic  energy  (KE)  rounds,  segmented  rods, 
multiple  launch  rocket  systems  (MLRS),  parachute  clusters,  submunition  dispersal,  and 
maneuvering  projectiles.  The  coupling  of  the  fluid  dynamic  solution  and  rigid-body  motion  is  a 
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major  advance  and  a  significant  accomplishment  that  eliminates  the  need  for  simplifying 
assumptions  and  provides  a  more  accurate  physically  based  simulation.  These  physically  based 
simulations  reveal  aerodynamic  flow-field  details  not  available  from  experiments,  screen  new 
concepts/designs,  and  minimize  expensive  flight  tests.  It  is  hoped  that  this  work  and  others  like 
it  will  create  a  change  from  the  test-simulate-test  environment  to  the  simulate-test-simulate 
environment. 
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